Nanoscale layering of antiferromagnetic and superconducting phases in Rb2Fe4Se5 
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We studied phase separation in a single-crystalline antiferromagnetic superconductor Rb2Fe4Se5 (RFS) us- 
ing a combination of scattering-type scanning near-field optical microscopy (s-SNOM) and low-energy muon 
spin rotation (LE-/iSR). We demonstrate that the antiferromagnetic and superconducting phases segregate into 
nanometer-thick layers perpendicular to the iron-selenide planes, while the characteristic in-plane size of the 
metallic domains reaches 10 flm. By means of LE-/iSR we further show that in a 40-nm thick surface layer the 
ordered antiferromagnetic moment is drastically reduced, while the volume fraction of the paramagnetic phase 
is significantly enhanced over its bulk value. Self-organization into a quasiregular heterostructure indicates an 
intimate connection between the modulated superconducting and antiferromagnetic phases. 

PACS numbers: 74.78.Fk, 76.75.+i, 68.37.Uv, 74.70.Xa 



The recent discovery of intercalated iron-selenide super- 
conductors has stirred up the condensed-matter com- 
munity accustomed to the proximity of the superconducting 
and magnetic phases in various cuprate and pnictide super- 
conductors. Never before has a superconducting state with a 
transition temperature as high as 30 K been found to coex- 
ist with such an exceptionally strong antiferromagnetism with 
Neel temperatures up to 550 K as in this new family of iron- 
selenide mateiials. The very large magnetic moment of 3.3 /Ib 
on the iron sites |6|, however, renders a microscopically ho- 
mogeneous coexistence of superconductivity and magnetism 
unlikely. Indeed, significant experimental evidence suggests 
that the superconducting and antiferromagnetic phases are 
spatially separated ||7 HT3]| . 

In the prototypical iron arsenide superconductors 
Ba(K)Fe(Co)2As2 both the phase separation lHH [151 and 
coexistence |T6l of antiferromagnetism and superconductivity 
have been shown to occur in certain regions of the phase dia- 
gram. At the same time, the structurally similar intercalated 
iron-selenide compounds (K,Rb,Cs)o gFei 6Se2 have defied 
all the efforts to synthesize a bulk single-phase material of this 
family. Absence of such electronically homogeneous super- 
conducting single crystals and a strong correlation between 
the superconducting and antiferromagnetic phases ifTTl [TSl 
necessitate a detailed research into the nature of their coexis- 
tence. The volume fraction of the magnetic phase has been 
estimated to 88% in recent Mossbauer [[TTl and bulk ^uSR 1 19 1 
studies. The shape of the phase domains, on the contrary, 
has seen much conflicting evidence with indications ranging 
from needlelike rather regular stripes |9| to insulating islands 
on a superconducting surface |8| to nanoseparated vacancy- 
disordered presumably metallic sheets in the bulk with an 
unknown in-plane form factor EOl . In a recent STM study 
of [110] KvFe2-ySe2 (KFS) thin films the superconducting 
phase was assigned to stoichiometric KFe2Se2 without iron 



vacancies f\2\. However, consistent theoretical description 
of the recent inelastic neutron scattering and angle-resolved 
photoemission measurements based on this assumption 
requires significantly different levels of the chemical potential 
in the bulk and at the surface (equivalent to a disparity of 
^ 0.08 electrons/Fe in the electron doping) i2n[22l . 

Infrared spectroscopy is well-suited to clearly distinguish 
between the antiferromagnetic semiconducting and paramag- 
netic metallic phases of RFS due to a large contrast in its 
complex dielectric function in this spectral range |TOl. To 
determine the geometry of the domains and unambiguously 
assign them to the semiconducting and metallic phases, a 
submicrometer-resolution technique must be used. In this Let- 
ter we employed apertureless, scattering-type scanning near- 
field optical microscopy (s-SNOM) in the infrared Il23lj25l . 
which enables determination of the material's complex di- 
electric function with an unsurpassed in-plane resolution of 
ca 20 nm and a typical nanometer topographic sensitivity of 
an atomic-force microscope (AFM). We complemented this 
spatially-resolved technique with low-energy muon spin rota- 
tion (LE-/iSR) Il26ll27l measurements on the same RFS single 
crystals to quantify the fraction of the magnetically-ordered 
phase in the bulk and trace its modification towards the sam- 
ple surface. 

The commercial s-SNOM near-field microscope uses an 
illuminated AFM probing tip to scan the sample surface 
and pseudo-heterodyne interferometric detection |28 1 to ex- 
tract the near-field amplitude and phase from the light scat- 
tered back from the tip. Standard platinum-coated AFM 
tips (Nano World ARROW-NCPt with a 25 nm radius) were 
used as antennas for nanofocusing of CO2 laser radiation at 
10.7 /zm (116 meV) wavelength (photon energy). The mea- 
surements were caiiied out on optimally-doped superconduct- 
ing RFS single crystals (batch BR26 in Ref.^ Tq « 32 K), 
cleaved prior to every scan. The sample surface was first char- 
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FIG. 1: (a) Microscope image of a 60 x 60 /xm^ surface patch of a fresiily-cleaved superconducting RFS single crystal. Typical rectangular 
15x8 flm^ area studied via near-field microscopy, (b) Superposition of the topography of a 15 x 8 /xm^ rectangular area (terrain) and the 
optical signal (brightness) normalized to that of silicon. Glossy areas indicate high silicon-RFS contrast and thus metallicity, while the matt 
areas are insulating. This combined response is broken down in (c) and (d) for the cross-section defined by the blue semi-transparent plane, 
(c) Optical contrast 52/5'2,si of the second harmonic ^2 of the near-field signal obtained at the 10.7 /im emission wavelength (116 meV photon 
energy) of a CO2 laser. Peaks in the contrast indicate metallic response (see text), (d) Displacement of the AFM tip while scanning along 
a 15 flm line obtained simultaneously with c. Dashed lines and blue shaded areas in (c) and (d) show the correlation between the metallic 
response and changes in the topography. 



acterized with a polarizing microscope. Figure [TJa) shows 
a 60 X 60 jxra? patch of the sample surface. A network of 
bright stripes is always observed on a freshly cleaved sur- 
face and does not depend on the polarization of the probing 
light. The stripes always occur at 45° with respect to the 
in-plane crystallographic axes. The topography of a repre- 
sentative 15 X 8 jim^ surface patch studied with an s-SNOM 
microscope is shown in Fig. [TJb) as a three-dimensional ter- 
rain image. One can clearly correlate the features in the AFM 
map with the bright stripes in the polarizing microscope im- 
age. This implies that the surface chemistry of this compound 
leads to inherent surface termination with mesoscopic terrac- 
ing upon cleaving, typically 10 — 30 nm high. Figures [Tfc) 
and [TJd) show the 2nd-harmonic near-field optical contrast 
(OC), obtained by normalizing the signal from the RFS sur- 
face to that that from a reference silicon surface (52/52, si), 
and the topography profile of RFS for the cross-section in- 
dicated with a blue translucent plane in Fig. [TJb). The am- 
plitude and phase of both the topography and the optical sig- 
nal are obtained simultaneously during a scan. Every peak in 
the OC maps signals a metallic optical response. The abso- 
lute values of the complex dielectric function were obtained 
within the extended finite-dipole model |30| using OC maps. 
The dielectric response of the semiconducting phase [dark re- 
gions in Fig.[TJb) and low OC in Fig.[TJc)] obtained in such a 
treatment of the experimental data (ei w 10, £2 ~ 0) is fully 
consistent with that of the single-phase semiconducting RFS 
crystals 1 10 |. The bright regions of the sample surface display 
negative values of ei , which provides solid evidence for their 
metallic character. 

By correlating the OC peak positions with the topography 
maps one can deduce that the location of the metallic regions 
on the sample surface is bound to the slopes in the terrain, as 



indicated in Figs.[T|c) and[T]J with dashed lines. Flat regions 
of the sample surface exhibit no metallic regions. Figure [TJd) 
shows that from the location of the OC peaks on the slopes 
one can identify well-defined sheets (horizontal dashed lines) 
approximately 10 nm apart. Metallic response (OC peak) is 
recorded whenever one of these sheets is exposed from the 
sample bulk at the intersections of the dashed lines with the 
topography line. Such correlation exists in all studied cross- 
sections of the OC and height maps. The thickness of the 
metallic sheets can be estimated by projecting the full width 
at half maximum (FWHM) of the OC peaks onto the surface 
topography [blue shaded areas in Fig. [TJc) and [TJd)] and on 
average amounts to 5 nm. The analysis of the coiTelations be- 
tween the OC peaks and the slopes in the sample topography 
carried out on multiple cross-sections of different near-field 
maps shows that the metallic sheets have an in-plane dimen- 
sion as well, which can be roughly estimated to w 10 jim in 
both directions. Therefore, phase separation in superconduct- 
ing RFS single crystals occurs on the nanoscale in the out-of- 
plane and mesoscale in the in-plane direction. From the thick- 
ness and periodicity of the metallic sheets we can estimate the 
volume fraction of the metallic domains in the 30 nm surface 
layer to approximately 50%, significantly larger than all pre- 
viously reported bulk values. 

It must be noted that the metallic volume fraction obtained 
by means of s-SNOM imaging only considers the near sur- 
face layer within « 30 nm. Whether or not this fraction 
changes with the distance from the sample surface and what 
the magnetic and superconducting properties of the metallic 
and semiconducting phases are cannot be determined from 
these measurements. Therefore, a microscopic sensor of the 
local magnetic moment with an adjustable implantation depth 
must be invoked. Such characteristics are provided by the 
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FIG. 2: (a) Schematic of the experiment with low-energy muons. The sample is depicted as a layered structure according to the discussion 
in the text, (b)-(d) Time dependence of /i+ spin polarization Af)P{t) in a zero magnetic field at 5 K after implantation in RFS single crystals 
at various depths determined by the muon stopping profiles of respective colors in (a), (e)-(g) Same for normalized /x+ spin polarization 
/para^(0 in a transverse magnetic field = 100 G, where the normalized asymmetry /p^fa is the paramagnetic volume fraction, (h)-(j) Long- 
time dependence of normalized fl^ spin polarization /para^'(0 in a transverse magnetic field //x = 100 G for the green muon stopping profile 
in (a) in the normal state at 46 K (h) and 32 K (i) and superconducting state at 5 K (j). Gray dashed lines in (e)-(j) show the slow relaxation 
envelope of the spin polarization; solid lines are fits to the data. 



LE-/iSR technique, which utilizes the predominantly spin- 
oriented positron decay of /i+ to detect the orientation of the 
muon spin. This orientation is influenced by the local mag- 
netic field at the muon implantation site, as well as magnetic 
fields applied externally. The experiment is shown schemati- 
cally in Fig. [2]; a). Based on the results of the near-field mea- 
surements the sample is sketched as an idealized periodic lay- 
ered structure, the finite in-plane dimension of the paramag- 
netic domains and the irregularity of their periodicity having 
been disregarded at this point. 

In the LE-juSR technique, in contrast to its conventional 
counterpart, the incident energetic MeV) muon beam is 
first moderated to about 15 eV in a condensed layer of solid 
N2 or Ar and then accelerated by a controlled electric field 
to achieve an adjustable implantation depth according to the 
muon stopping profile in a given material, calculated using 
the Monte Carlo algorithm TRIM.SP \3T\. To study the de- 
pendence of the phase separation in RFS single crystals on 
the distance from the sample surface, we chose three muon 
stopping profiles shown as a red (achieved by accelerating the 
moderated muon beam to the energy of 1 keV), green (9 keV), 
and blue (14 keV) shaded profile in Fig. |2ja). The resulting 
time dependence of jJ.^ spin polarization AoP(f) at zero exter- 
nal magnetic field is shown in Figs.|2|b)-(d), the colors coiTe- 
spond to those of the muon stopping profiles in Fig.|2|a). As is 
clear from Fig.|2|a), the shallowest muon stopping profile (red 
shaded area) is peaked at about 10 nm and extends to about 
30 nm below the sample surface. It thus probes a region of the 
sample bulk comparable to that observed via s-SNOM. The 
relaxation of the muon polarization after implantation in the 
sample occurs at two different time scales: fast depolarization 



takes place during the first 150 — 250 ns and is associated with 
the ordered antiferromagnetic phase, while a much slower 
evolution dominates thereafter and stems predominantly from 
the paramagnetic volume fraction. The fast depolarization 
rate is proportional to the width of the magnetic-field dis- 
tribution at the muon stopping site in the ordered phase and 
thus to the antiferromagnetic moment. Figures |2|c) and |2|d) 
show that in the bulk [green and blue muon stopping profiles 
in Fig. |2|a)] the fast component is approximately the same 
and rather narrow, whereas near the surface it broadens sig- 
nificantly [Fig. [2|b)]. Using a two-component fit model for 
the zero-field /iSR data, the fast depolarization rates in the 
bulk were found to agree within the error bars: 52(13) /zs and 
37(7) /is for the green and blue profiles, respectively, whereas 
it is significantly reduced to 21 (4) /is close to the surface (red 
profile). Under the assumption that the muon stopping sites 
in the unit cell are the same for all three profiles and that the 
depolarization rate should be the same in the bulk, the antifer- 
romagnetic moment is reduced to only 50% of its bulk value 
in the 30 nm surface layer. 

The muon decay asymmetry in Figs.[2|b)-(d) indicates that 
the paramagnetic volume fraction (slow component) is en- 
hanced closer to the surface. To avoid the contaminating 
/iSR signal originating from muon decays in the Ni sample 
holder, we caiTied out transverse-field measurements in the 
same configuration and at the same temperature of 5 K as 
in Figs. |2|b)-(d). In this case the oscillating component of 
the muon decay asymmetry comes only from the paramag- 
netic domains in the sample and can be normalized accord- 
ingly to become fparaP{t), where /para is the paramagnetic 
volume fraction, given by the value of the slowly depolariz- 
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ing envelope of the muon spin precession (gray dashed lines 
in Fig. [2|i at zero time. From Figs, ^f) and|2|g) one imme- 
diately infers that well in the sample bulk the paramagnetic 
phase constitutes about 20% of the sample volume. It thus 
characterizes the phase separation in the bulk of the supercon- 
ducting RFS single crystals as probed by other techniques, 
such as Mossbauer spectroscopy and conventional /zSR. On 
the other hand, Fig.|2|e) shows that in the 30 nm surface layer 
[red stopping profile in Fig. [2|a)] the paramagnetic volume 
fraction strongly increases to ss 50%. Such a high value most 
likely occurs due to a reduced constraining potential between 
antiferromagnetic sheets with a significantly smaller ordered 
moment close to the sample surface, which leads to an expan- 
sion of the metallic regions. These depth-dependent physical 
properties could explain the difference in the chemical poten- 
tial required for a consistent interpretation of the data obtained 
with bulk (inelastic neutron scattering) and surface-sensitive 
(ARPES) probes. 

To characterize the superconducting phase of RFS with 
/iSR we studied the temperature dependence of the oscillat- 
ing component of the muon asymmetry in a vortex state gen- 
erated by the transverse external magnetic field, perpendicu- 
lar to the metallic sheets. The oscillations are damped by the 
inhomogeneity of the magnetic field in the vortices, which, in 
turn, is proportional to the condensate density ris °^ l/'^ab- 
measurements were carried out in a transverse magnetic field 
H± = 100 G for the muon stopping profile peaked at 80 nm 
below the sample surface [green shaded area in Fig. |2|a)]. 
The time dependence of the normalized /z^ spin polariza- 
tion /para^(0 in the normal state at 46 and 32 K and in the 
superconducting state at 5 K is shown in Figs. |2]^h)-(j), re- 
spectively. It is evident from the data that the damping in 
the normal state is approximately constant, while in the su- 
perconducting state it is noticeably faster, which indicates the 
presence of a superconducting condensate. By subtracting the 
normal-state damping from that in the superconducting state 
we can estimate the London penetration depth Ay w 550 nm. 
Assuming that the phase separation in the bulk has the same 
layered structure as detected near the surface, this value must 
be revisited in an appropriate model. For a stack of supercon- 
ducting layers of thickness dsc separated by insulating layers 
of thickness d^^s the inhomogeneity of the magnetic field in a 
vortex (which is then a stack of two-dimensional vortex pan- 
cakes) is described by the Lawrence-Doniach in-plane pene- 
tration depth related to the bulk London penetration depth 

A,ab via A|| ~ Xab{dsc +<^ins)'''^/<^s/^ (see Ref.|32j). This reduc- 
tion reflects the fact that A| | is related to the average supercon- 
ducting condensate density (nj) = n^d^^/ (c/sc +^^ins) = «s/paia- 
Taking these considerations into account and using the bulk 
paramagnetic volume fraction /para ~ 0.2 obtained in our LE- 
/zSR measurements one can estimate the intrinsic bulk in- 
plane London penetration depth of the superconducting phase 
to be Aab ~ 250 nm. This value would increase towards A| | due 
to the finite in-plane dimension of the superconducting layers 
and their disordered stacking since the superconducting phase 



inclusions would then contribute independently to the /xSR 
depolarization. At the same time, our procedure slightly over- 
estimates the value of A|| due to the widening of the vortex 
field distribution close to the surface ||33l . 

Similar effective-medium approximations (EMA) must be 
used for an adequate analysis of the results obtained with other 
experimental techniques. A recent study of the optical con- 
ductivity of the same superconducting single crystals ifTOl em- 
phasized the importance of EMA but could not make a con- 
crete estimate due to the unknown details of the phase separa- 
tion in this compound. It reported the total plasma frequency 
of the itinerant charge carriers of about 100 meV. Using the 
same dimensions for the metallic paramagnetic and semicon- 
ducting antiferromagnetic layers as in the interpretation of the 
LE-/ZSR results one can extract the bulk superconducting op- 
tical response from the conductivity data reported in Ref.fTOl 
A simple fit of the experimental data as the optical response 
of a perfect superlattice gives a significantly larger value of 
the total plasma frequency OJpj' w 300 meV and brings the es- 
timated value of the London penetration depth A°,^' ~ 2 /zm 
closer to that obtained with LE-/iSR but still 8 times larger 
Accounting for the finite in-plane dimension of the paramag- 
netic domains is expected to lead to a better agreement be- 
tween the two techniques. 

The origin of the phase separation observed in this work can 
lie either in a chemical stratification into e.g. iron-vacancy 
ordered (antiferromagnetic) and disordered (paramagnetic) 
phases or in a purely electronic segregation on a homoge- 
neous crystalline background. Self-organization of a chem- 
ically homogeneous structure into (quasi)periodically segre- 
gated phases is not unprecedented — a similar phenomenon 
has been observed in copper-based superconductors, where 
antiferromagnetic stripes of copper spins were found to be 
spatially separated by periodic domain walls close to a partic- 
ular hole doping level of 1/8 |34|. Be it of chemical or elec- 
tronic nature, the phase separation in superconducting RFS 
single crystals reported here represents an interesting case of 
a naturally-occurring quasi-heterostructure. 

To conclude, by combining the unique optical imaging ca- 
pabilities and nanoscale resolution of the s-SNOM near-field 
microscope with bulk sensitivity at variable depth of LE-/xSR 
we determined the geometry and magnitude of the phase sep- 
aration in RFS superconducting single crystals. The param- 
agnetic domains were found to have a shape of thin metallic 
sheets parallel to the iron-selenide plane of the crystal with 
a characteristic size of only several nanometers out of plane 
but up to 10 pLm in plane. By means of LE-/xSR we further 
show that the antiferromagnetic semiconducting phase occu- 
pies w 80% of the sample volume in the bulk and is strongly 
weakened near the surface. These results have important im- 
plications for the interpretation of bulk- and surface-sensitive 
measurements on Rb2Fe4Se5, and for the understanding of the 
interplay between superconductivity and antiferromagnetism 
in this material. 
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